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Abstract 

Using a sample of 1.06 x 10 8 if/ decays collected by the BESIII detector, XcO and Xc2 decays 
into 7r°7r° and r/r/ are studied. The branching fraction results are Br(xc0 — > vr°7r°) = (3.23 ± 
0.03 ± 0.23 ± 0.14) x 10" 3 , Br( Xc 2 -> vr 7r°) = (8.8 ± 0.2 ± 0.6 ± 0.4) x 10~ 4 , Br( Xc o -> w) = 
(3.44 ± 0.10 ± 0.24 ± 0.2) x 10 -3 , and 5r(x C 2 -> w) = (6.5 ± 0.4 ± 0.5 ± 0.3) x 10" 4 , where 
the uncertainties are statistical, systematic due to this measurement, and systematic due to the 
branching fractions of ip' — > 7X c j, respectively. The results provide information on the decay 
mechanism of Xc states into pseudoscalars. 

PACS numbers: 13.25.Gv, 14.40.Pq 
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I. INTRODUCTION 



In the quark model, the Xcj {J — 0, 1, 2) mesons are L = 1 cc states. Since they cannot 
be produced directly in e + e~ collisions, they are not as well studied as the ip states. On the 
other hand, ip' — > ^Xcj decays yield many Xcj mesons, providing a clean environment for 
Xcj investigations. In this paper, we study two-body decays of the Xco and Xc2 into 7r°7r° and 
rji] final states 1 . Knowledge gained from these decays provides information on both the Xcj 
parents and their pseudo-scalar daughters, as well as a greater understanding of the decay 
mechanisms of Xcj mesons [1]. 

Recently, v^n and Xc2 decays into two-meson final states were studied by the CLEOc 
collaboration [2|. In this analysis, we use a sample of 1.06 x 10 8 ip' decays collected by the 
BESIII detector to perform a study of these decays. 



II. BESIII AND BEPCII 

The analysis reported here is based on about 1.06 x 10 8 ip' events collected by BESIII 
at BEPCII. BEPCII/BESIII is a major upgrade of the BESII experiment at the BEPC 
accelerator J4] for studies of hadron spectroscopy and r-charm physics j^]. The design peak 
luminosity of the double-ring e + e~ collider, BEPCII, is 10 33 cm" 2 ^" 1 at a beam current of 
0.93 A. The BESIII detector with a geometrical acceptance of 93% of Att, consists of the 
following main components: 1) a small-celled, helium-based main draft chamber (MDC) with 
43 layers. The average single wire resolution is 135 /im, and the momentum resolution for 
1 GeV/c charged particles in a 1 T magnetic field is 0.5%; 2) an electromagnetic calorimeter 
(EMC) made of 6240 Csl (Tl) crystals arranged in a cylindrical shape (barrel) plus two 
endcaps. For 1.0 GeV photons, the energy resolution is 2.5% in the barrel and 5% in the 
endcaps, and the position resolution is 6 mm in the barrel and 9 mm in the endcaps; 3) a 
Time-Of-Flight system (TOF) for particle identification composed of a barrel part made of 
two layers with 88 pieces of 5 cm thick, 2.4 m long plastic scintillators in each layer, and 
two endcaps with 96 fan-shaped, 5 cm thick, plastic scintillators in each endcap. The time 
resolution is 80 ps in the barrel, and 110 ps in the endcaps, corresponding to better than 
a 2 sigma K/pi separation for momenta below about 1 GeV/c; 4) a muon chamber system 



1 We do not consider \ci decays into these final states, as they are forbidden by spin-parity conservation. 
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(MUC) made of 1000 m 2 of Resistive Plate Chambers (RPC) arranged in 9 layers in the 
barrel and 8 layers in the endcaps and incorporated in the return iron of the superconducting 
magnet. The position resolution is about 2 cm. 

The optimization of the event selection and the estimation of physics backgrounds are 
performed through Monte Carlo simulations. The GEANT4-based simulation software 
BOOST [G] includes the geometric and material description of the BESIII detectors, the 
detector response and digitization models, as well as the tracking of the detector running 
conditions and performance. The production of the if)' resonance is simulated by the Monte 
Carlo event generator KKMC [7], while the decays are generated by EvtGen [8[ for known 
decay modes with branching ratios being set to the PDG 10] world average values, and by 
Lundcharm 

[9] 

for the remaining unknown decays. The analysis is performed in the frame- 
work of the BESIII Offline Software System (BOSS) [ll|] which takes care of the detector 
calibration, event reconstruction and data storage. 



III. EVENT SELECTION 



A photon candidate is defined as a shower in the EMC with an energy deposit exceeding 
50 MeV. The tt° and rj candidates are reconstructed from pairs of photon candidates, using 
the average event vertex of each run as the assumed origin of the photons. For 7r° — > 77, 
the 77 invariant mass is required to satisfy 0.075 GeV/c 2 < M(jj) < 0.175 GeV/c 2 . For 
rj —> 77, the 77 invariant mass is required to satisfy 0.458 GeV/c 2 < M(77) < 0.608 GeV/c 2 . 
The decay angle of a photon is the polar angle measured in the ir° or 77 rest frame with respect 
to the 7T° or rj direction in the if)' rest frame. Real 7r° and 77 mesons decay isotropically, and 
their angular distributions are flat. However, the 7r° and 77 candidates that originate from a 
wrong photon combination do not have a flat distribution in this variable. To remove wrong 
photon combinations, the decay angle is required to satisfy | cos 9d eC ay\ < 0.95. 

Candidate events for the final states of interest (77r°7r° and 77777 ) are selected using 
the following basic selection criteria. An event must have 5 or 6 photons and no charged 
tracks. All possible two photon pairings (the radiative photon from the if>' decay which has 
E < 0.4 GeV is excluded) in the event are used to form 7r° and 77 candidates. The candidate 
event uses the photon pairings giving the minimum 
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with Pi and P2 being the pulls, defined as: 



P(tt%) 



M, 



77 



I* 



Hi 



where M 77 is the reconstructed 77 invariant mass, m n o/ v is the known 7r° or 77 mass IOJ, and 



o~ 77 is the 77 mass resolution, with typical values of 7 MeV/c 2 for the 7r° and 12 MeV/c 2 
for the 77. If there is more than one radiative photon candidate (E<0.4 GeV), the one that 
gives the least |M 57 — is used. 

Backgrounds with missing final state particles are suppressed by requiring small trans- 



verse momentum squared p\ 



; 7 J 



Pt-y = ^Pmiss Sm 2 (# 7 /2), 



where p m i SS is the missing momentum opposite to the n n or r\r\ system and 9 1 is the angle 
between the radiative photon and the direction of the missing momentum p m iss- The 77r 7r° 
events are required to satisfy p 2 7 < 0.04 (GeV/c) 2 , while the 77777 events are required to 
satisfy p 2 tl < 0.01 (GeV/c) 2 and x V v < 4 - 

To study the efficiency of the if)' — > •JXcj, Xcj 7r°vr° and Xcj VV selection, Monte 
Carlo samples for each Xcj state into each final state are generated using a (1 + Acos 2 #) 
distribution, where 9 is the radiative photon angle relative to the positron beam direction, 
and A = 1 for Xco an d A = 1/13 for x c2 , in accordance with expectations for El transitions. 
The decay products of the Xco are generated using a flat angular distribution, while those of 
the Xc2 are generated according to a double correlation function of the polar angles of the 



mesons measured in the Xc rest frame relative to the transition photon direction 
The efficiencies obtained from the Monte Carlo simulation are shown in Table [H 



12| 



TABLE I: Efficiencies (in %) obtained from analysis of Monte Carlo generated events. 



Mode 


XcO 


Xc2 


7r°vr 
rjf] 


55.6 ±0.2 
40.3 ±0.2 


59.8 ±0.2 

43.9 ±0.2 



Figured] and Figure [2] show comparisons in the Xco region between data and Monte Carlo 
simulation for the selection criteria used. The good agreement across the distributions shows 
that the efficiency estimated from Monte Carlo simulation is reliable. 
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FIG. 1: Comparisons between data and Monte Carlo simulation of ip' — > ^Xcj,XcJ 7r°7r° for 
selection criteria used, (a) Photon multiplicity distribution, (b) The 77 invariant mass distribution 
for 7T° candidates, (c) The distribution ofp^. (d) The X-k°tt° distribution. Dots with error bars are 
data in the XcO region. The histogram is the Monte Carlo simulation for ip' — > JXc0,Xc0 7r°7r° 
plus normalized background estimated from inclusive ip' Monte Carlo samples. 



IV. BACKGROUND ANALYSIS 



The backgrounds in the selected event sample from a number of potential background 
channels listed in the PDG [lO] are studied with Monte Carlo simulations. The main back- 
ground to Xcj — > vr 7r° originates from ip' — > ^Xcj, Xcj — > iJ/'fp, J/ip — > IV- Using the world 
average branching fractions [lO] for this mode, we estimate that 48 events from this channel 
are in the signal region. However, the simulation also shows that the background does not 
peak at the Xco nor the Xc2 mass region. The main backgrounds to Xcj VV originate from 
ip' — > ir 7r°J/ip and ip' — > r)J/ip, J/ip — > 777. There are about 233 surviving background 
events in the signal region. 




0.02 0.04 0.06 0.08 0.1 

pf 7 (GeV 2 /c 2 ) 
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FIG. 2: Comparisons between data and Monte Carlo simulation for ip' —> jXcJ,XcJ VV f° r the 
selection criteria used, (a) Photon multiplicity distribution, (b) The 77 invariant mass distribution 
for r) candidates, (c) The distribution of pL. (d) The Xrm distribution. Dots with error bars are 
data in the XcO region. The histogram is the Monte Carlo simulation for ip' — >■ 7Xc0i XcO — > f]f] plus 
the normalized background estimated from inclusive ip' Monte Carlo samples. 

A 10 8 inclusive ip' Monte Carlo event sample is also used to investigate other possible 
surviving background events. Figs. [3] (a) and [3] (b) show the radiative photon energy dis- 
tribution of the selected Xcj — > vr°7r° and \ c j — > rjr] events and the normalized backgrounds 
estimated with the inclusive ip' Monte Carlo sample, respectively. In the Xcj signal region, 
there is no peaking background from the inclusive ip' Monte Carlo sample. 

The background in our signal region originating from non-resonant processes is studied 
using a continuum data sample collected at a center of mass energy of 3.65 GeV. Normalized 
according to the luminosities, the contribution to Xcj ~ > 7r°vr° is 384 events, as shown in 
Fig. [3] (a), and the contribution to Xcj — > VV is 48 events, as shown in Fig. |3] (b). These 
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FIG. 3: Radiative photon energy distributions of (a) selected Xc 7r 7r° events, and (b) selected 
Xc —> VV events. Dots with error bars are data. The open histogram is the normalized back- 
ground estimated from the inclusive ip' Monte Carlo sample and from the continuum. The shaded 
histogram is the normalized contribution from the continuum. 



backgrounds are small, do not peak in the signal region, and are included as part of the 
polynomial background below. 



V. NUMBER OF ip' EVENTS 

The number of if)' events, N$>, used in this analysis is determined from the number of 
inclusive hadronic if)' decays. Charged tracks are selected requiring their point of closest 
approach to the beam axis be within 1 cm of the beam line, and their angle with respect to 
the beam axis, 9, must satisfy | cos#| < 0.93. Photon candidates must have at least 25 (50) 
MeV of energy in the barrel (end-cap) EMC, and have | cos#| < 0.93. 

Event selection requires at least one charged track. To remove beam associated back- 
ground and background from Bhabha events, there are special requirements on low charged 
multiplicity events. For events with one charged track, there must be at least three photons, 
the acolinearity angle between the two highest energy photons must be greater than 7°, and 
the total energy in the EMC, E EMC , must be greater than 0.2 of the center of mass energy, 
E cm , and less than 0.85E cm . Events with two or three tracks must have Eemc > 0.2E cm 
in order to suppress beam associated backgrounds. Backgrounds from Bhabha events are 
reduced by requiring the presence of at least two photons and Eemc < 0.85E cm or the 
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largest energy deposit in the calorimeter less than 0.85 times the beam energy, E^ eam . In 
addition, the second highest momentum track must have momentum less than 0.9Eb ea m, 
and the acolinearity angle in the x — y plane of the two highest momentum tracks must be 
greater than 7°. 

The number of hadronic events is determined from the distribution of z, which is the 
average of the distances, z, from the interaction point along the beam of the point of clos- 
est approach of tracks to the beam line. Two methods are used: fitting the distribution 
with a Gaussian plus a 2nd order polynomial background and counting events in a signal 
region and subtracting sideband events. Backgrounds from Bhabha, dimuon, and ditau 
events surviving the selection criteria are very small. The continuum contribution and the 
surviving backgrounds are removed by subtracting the number of events selected with the 
above criteria from a continuum sample taken at a center of mass energy of 3.65 GeV and 
normalized by relative luminosity and the 1/s dependence. The efficiency for ip' — > hadrons 
is determined by simulation p| and is 0.80. The agreement between data and Monte Carlo 
(MC) simulation is shown for the distribution of the number of charged tracks in Fig. H] (a) 
and for Eemc in Fig- HI (b). 




5 10 15 0.5 1 1.5 



N trk Total E^ MC /E cm 

FIG. 4: (a) The distribution of the number of charged tracks for events satisfying selection cri- 
teria, (b) The distribution of the total energy in the EMC divided by the center of mass energy, 
Eemc / E C m, f° r events satisfying selection criteria. All requirements are applied to events with one 
to three charged tracks except the EMC requirements. Dots are data, the light shaded histogram 
is the sum of normalized continuum and V' hadrons simulated events, and the dark shaded 
histogram is from continuum data. 
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The result is N^> = (1.06 ± 0.04) x 10 8 , where the error is systematic and is determined 
mostly by the track efficiency difference between data and MC (1.2%), the variation with 
the minimum charged track multiplicity requirement (2.86%), the difference when a mini- 
mum transverse momemtum requirement is used (0.95%), the uncertainty of the generator 
model (0.61%), and error due to the continuum subtraction (0.91%). The statistical error is 
negligible. A second analysis using a much different selection criteria with a higher efficiency 
determines an almost identical result. 



VI. FITTING RESULTS 

The Xc —> vr 7T branching fraction is calculated using 



0\ Nobs 



Njf,/ ■ e ■ Br(ip' — > 7Xcj) • Br(ir° — > 77) • Br(ir° — > 77) ' 

where N b s is the number of events observed, N^> is the number of ip' events, and e is the 
selection efficiency obtained from MC simulation. 

The radiative photon energy spectrum of Xcj vr°7r° candidates, shown in Fig. [5J is 
fitted using an unbinned maximum likelihood fit in the range from 0.06 GeV to 0.36 GeV. 
The shapes of the Xco an d Xc2 are obtained from Monte Carlo simulation and the masses 



and widths of Xcj are fixed to their PDG values [lOj. A 2nd-order Chebychev polynomial is 
used to describe the backgrounds, including those found in the inclusive MC study and the 
continuum. The fit gives a Xco signal yield of 17443 ± 167 events and a Xc2 signal yield of 
4516 ± 80 events. The selection efficiency from Monte Carlo simulation of ip' — > 7Xco(Xco 
7T 7r°,7r — > 77) is (55.6 ± 0.2)% and the efficiency of t])' —> r fXc2(Xc2 7r°7r ,7r° — > 77) is 
(59.8 ± 0.2)%. The branching fractions are then determined to be 

Br(Xco -> 7T 7r°) = (3.23 ± 0.03) x 10" 3 , 

Br(xc2 ->■ vr°7r°) = (8.8 ± 0.2) x 10" 4 , 

where the errors are statistical only. 

The fit to the radiative photon energy spectrum of Xcj ~^ VV candidates, shown in Fig. |6l 
gives a Xco signal yield of 2132 ± 60 events and a Xc2 signal yield of 386 ± 25 events. The 
selection efficiency is 40.3 ± 0.2% and 43.9 ± 0.2% for Xco VV an d Xc2 Wit respectively. 
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E(y) (GeV) 

FIG. 5: The radiative photon energy spectrum of selected \ c — > 7r 7r° events. Dots with error bars 
are data. The solid curve is the result of a fit described in the text. The dotted curves are the \cj 
signals. The dashed curve is the background polynomial. 




E(y) (GeV) 

FIG. 6: The radiative photon energy spectrum of selected \c W events. Dots with error bars 
are data. The solid curve is the result of a fit described in the text. The dotted curves are the Xcj 
signals. The dashed curve is the background polynomial. 

The branching fractions are 

Br{ X co -> VV) = (3-44 ± 0.10) x 10~ 3 , 

Br( X c2 -> m) = ( 6 - 5 ± °-4) x 10~ 4 , 
where the errors are statistical only. 
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VII. SYSTEMATIC UNCERTAINTIES 



The systematic uncertainties on the branching fractions come from many different sources 
and are summarized in Table [III The uncertainty due to photon detection and photon 
conversion is 1% per photon. This is determined from studies of photon detection efficiencies 
in well understood decays such as J /if) — > p°7r° and study of photon conversion via e + e _ — > 
77 • 

The uncertainty due to 7r° selection is determined from a high purity control sample of 
J/ip — > 7T + 7r~7r° decays. The ir° selection efficiency is obtained from the change in the 7r° 
yield in the 7r + 7r~ recoiling mass spectrum with or without the 7r° selection requirement. 
The difference of it reconstruction efficiency between data and Monte Carlo simulation 
gives an uncertainty of 1% per 7r°. The uncertainty from the rj selection is 1% per rj, which 
is determined in a similar way from a high purity control sample of J /if) — > rjpp. 

The systematic error from the p% requirement is determined by not using the requirement. 
The change in the yield gives systematic errors of 0.9% for Xco vr°7r°, 1.2% for Xc2 — > vr°7r°, 
0.1% for Xco Wit and 0.3% for Xc2 VV- The uncertainties from Xm requirement are 
0.6% for Xco ~^ VV and 2.6% for Xc2 VV) and are determined in a similar way. 

Since the shapes of the signals in the fit are obtained from Monte Carlo simulation, 



their uncertainties are estimatec 
deviation from the PDG values 



by changing the masses and widths of Xcj by one standard 



10j and taking into account the uncertainties of the photon 
energy scale and resolution in the Monte Carlo simulation. They are 1.6% for Xco — >■ vr 7r°, 
1.2% for Xc2 vr 7T°, 1.4% for Xco ~^ VV) an d 1-5% Xc2 — > VV- 

The background uncertainties are evaluated by changing the background fitting function 
from a second order polynomial to third order, resulting in changes of branching ratios by 
0.5% for Xco -> 7r 7T°, 0.5% for Xc2 -> 7r°7r°, 0.2% for Xco -> VV) and 0.3% for Xc2 -> VV- 

The systematic uncertainties due to the fitting of the radiative photon energy spectrum 
were evaluated by changing the fitting range from (0.05, 0.37) GeV to (0.07, 0.35) GeV. 
The change in yield for this variation gives systematic uncertainties of 0.3% for Xco 7r 7r°, 
0.3% for Xc2 vr°7r°, 0.8% for Xco ~^ VV) and 1.3% for Xc2 VV- 

The systematic uncertainties due to the trigger efficiency in these neutral channels is 
estimated to be < 0.1%, based on cross checks using different trigger conditions. The 
uncertainty on the number of if)' events is 4%. 
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The total systematic uncertainties, shown in Table. [TTJ are obtained by adding all the 
above systematic errors in quadrature. The uncertainty due to the ip' — > 7% c branching 
fractions is kept separate and quoted as a second systematic uncertainty. 



TABLE II: Systematic uncertainties expressed in percent. 



Mode 


Y^.n — > 7T 7r 
AcU ' " n 


Y„o — > 7T 7r° 
Acz ' " " 


Y^.n — > 7777 
AcU ' 77 


AcZ 7 77 


photon detection 


5 


5 


5 


5 


7t®(ti) reconstruction 


2 


2 


2 


2 


7^17 


0.9 


1.2 


0.1 


0.3 


Xrp] 






0.6 


2.6 


signal shape 


1.6 


1.2 


1.4 


1.5 


background shape 


0.5 


0.5 


0.2 


0.3 


fitting range 


0.3 


0.3 


0.8 


1.3 


trigger 


0.1 


0.1 


0.1 


0.1 


AV 


4 


4 


4 


4 


Total 


7.0 


6.9 


6.9 


7.5 



TABLE III: Branching fraction results (in units of 10 -3 ) for each decay mode. The uncertainties 
are statistical, systematic due to this measurement, and systematic due to the branching fractions 
of ip' — > 7Xcjj respectively. CLEOc results are determined using their own branching fractions 
for ijj' —7- 7Xcj; while ours are determined using branching fractions from the PDG. If we use the 
CLEOc branching fractions, we find Br(xc0 —> 7r°7r°) = 3.29 x 10~ 3 , Br(xcO -> w) = 3-51 x 10" 3 , 
Br(xc2 -> vr°7T ) = 0.78 x 10~ 3 , and Br(x C 2 -> w) = °- 58 x 1°~ 3 - 



Mode 




XcO 


Xc2 


TT 7T 


This Work 


3.23 ±0.03 ±0.23 ±0.14 


0.88 ±0.02 ±0.06 ±0.04 




CLEOc [2] 


2.94 ±0.07 ±0.32 ±0.15 


0.68 ±0.03 ±0.07 ±0.04 




PDG [10J 


2.43 ±0.20 


0.71 ±0.08 


TjTj 


This Work 


3.44 ±0.10 ±0.24 ±0.13 


0.65 ±0.04 ±0.05 ±0.03 




CLEOc [2] 


3.18 ±0.13 ±0.31 ±0.16 


0.51 ±0.05 ±0.05 ±0.03 




PDG [10] 


2.4 ±0.4 


< 0.5 
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VIII. SUMMARY 



In summary, with a sample of 1.06 x 10 8 if)' events in the BESIII detector, improved 
measurements of the branching fractions of Xco,2 vr°7r° and Xco,2 ~* VV are performed: 
Br(Xco -> vr 7T°) = (3.23 ± 0.03 ± 0.23 ± 0.14) x lO" 3 , Br( X c2 -> vr°7r°) = (8.8 ± 0.2 ± 
0.6 ± 0.4) x 10~ 4 , Br(xco -> w) = (3-44 ± 0.10 ± 0.24 ± 0.2) x 10^ 3 , and Br( X c2 ->■ w) = 
(6.5 ± 0.4 ± 0.5 ± 0.3) x 10~ 4 , where the uncertainties are statistical, systematic due to this 
measurement, and systematic due to the branching fractions of if)' — > 'JXcj, respectively. 
Results are listed in Table II III and compared with previous measurements. 
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